Two-level quantum dynamics, integrability and unitary NOT gates 
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We study the dynamics of a two-level quantum system interacting with an external electromag- 
netic field periodic and quasiperiodic in time. The quantum evolution is described exactly by the 
classical equations of motion of a gyromagnet in a time-dependent magnetic field. We prove that 
this classical system is integrable as a consequence of the underlying unitary quantum dynamics. 
As a consequence, for the periodic case: i) rigorous assessment of the validity of the rotating-wave 
approximation (RWA) becomes possible even beyond the assumptions of resonance and weak cou- 
pling (the latter conditions are also shown to follow from the method of averaging); ii) we determine 
conditions for the realization of the quantum NOT operation beyond the RWA, by means of classi- 
cal stroboscopic maps. The results bear upon areas as diverse as quantum optics, nuclear magnetic 
resonance (NMR) and quantum computation. 
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Two-level systems have achieved a paradigmatic place 
in modern quantum mechanics with the advent of strong 
laser pulses and advanced NMR techniques. By coupling 
a two-level system to a time-dependent classical field it 
is possible to describe the equivalent dynamics of a two- 
level atom interacting with an external electromagnetic 
field (quasi)pcriodic in time in the dipolc approximation 
(relevant to quantum optics) or of a nuclear spin 1/2 in an 
external (quasi)pcriodic magnetic field (relevant to NMR 
and to the NOT operation in quantum computation). 
They are modelled by the same Hamiitonian 

H(t) = ~B(t) ■ E, (1) 

where S = (a x ,a y ,a z ) is the vector composed by the 
Pauli matrices. For simplicity we adopted H = 1 and 
thus the magnetic field is in units of frequency. We shall 
analyze two important external periodic fields frequently 
used in quantum optics and NMR experiments, particu- 
larly for the implementation of the quantum NOT gate. 
The first one is a radio-frequency field Q given by 

B R = -2(B cos (Jt,B sin (Jt,Bi), (2) 

where Bq and B\ are constant amplitudes. We call it 
rotating (R) because it corresponds to a magnetic field 
which rotates around the z-axis. In quantum optics it 
corresponds to a circularly polarized radiation field. The 
dynamics for this field may be solved analytically by ap- 
plying the unitary transformation exp(— iu>ta z ) to the 
wave-function in Schrodinger equation for the Hamiito- 
nian {Tjl: in the resulting rotating frame the Hamiitonian 
becomes time-independent and the solution emerges eas- 
ily 0] • We also investigate the dynamics under a periodic 
nonrotating (NR) magnetic field along a preferred direc- 



tion: 

B NR = -2(B 2 coscot,0,B 3 ), (3) 

with B2 and B3 time- independent amplitudes. Bjv/j cor- 
responds to a plane-polarized radiation field in quantum 
optics 0, a fact which explains its conceptual impor- 
tance. In fact, the quantum dynamics corresponding to 
BpfR is very rich and far from soluble: it has been a 
recurrent theme in quantum optics since the seminal pa- 
pers of Bloch and Siegert Q and Autlcr and Townes |]| . 
This fact has often led to the replacement of Bn R by 
Bj{, in particular for the dynamics of the quantum com- 
putation However, the R field is much more difficult 
to manufacture and, therefore, the most commonly used 
in NMR experiments, in particular for quantum compu- 
tational purposes, is Bnr Thus, in each of the above 
areas, quantum optics, NMR and quantum computation, 
it is of central importance to understand the dynamics 
induced by B^r. 

The above-mentioned replacement of B^r by B^ is 
called rotating-wave approximation (RWA). It has been 
heuristically discussed in many text books (see for the 
most careful presentation), but we do not know any ex- 
act statement on the domain of validity of the various 
parameters (in contrast, an exact method, which may 
be called rotating-wave transformation, has been intro- 
duced in 0). Such a statement may be obtained from 
the following result. 

Result 1 : Let B 2 T <C 1 (weak coupling), with T = 2-k/uj. 
Then, if u = 2B% (resonance), the solution of the 
Schrodinger equation tJjr corresponding to H(t) given by 
(|TJ) with B = B R (with B Q = B 2 /2 and B x = B 3 ) ap- 
proximates the solution ipNR of the Schrodinger equation 
with the Hamiitonian Q with B = Bjv/j in the sense that 
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||^fl(f)-^j\rJt(t)|| - 0(B 2 T), whenever < t < O^). 

Above, 1 1 • 1 1 is the Euclidean norm in the spinor vector 
space. Writing 



ip(t) = cxp 



iujta z 



(4) 



we obtain for ip' under the condition of resonance the 
Schrodingcr equation 



(5a) 



where e = B 2 to 1 is, by assumption, a small parameter, 
and 



— (a x cos 2ojt — (Ty sin 2u)t) ip 1 (5b) 



is periodic in t, of average equal to ^-^-tp'. By the aver- 
aging theorem (||, Theorem 11.1), the solution ip' of (JSJ 
is 0(B 2 T)-close to the solution of 



#6 
1 dt 



B 2 lo' 



(6) 



if < i < C(-g|y); undoing the transformation |0J, 
we obtain that "0 is close to the function ipn(t) = 
e -iwt CTz /2 e - l B 2 ia a> /2^ R ^Q^ wmc h j s t he wave-function in 

the RWA. 

A different approach has been pioneered by Feynman 
et al (see also and 0). Considering a density 
operator p satisfying the conditions of unitary evolution 
(ip = [H(t),p]), normalization (Trp = 1) and purity 
(p 2 = p; this condition is, however, not essential to the 
method), one may define the paramctrization 



S • S) 



(7) 



to map the quantum dynamics in terms of classical 
Hamiltonian formalism |9l. defined by 



H{t) = -B(f) • S, 
S 2 = l. 



(8a) 
(8b) 
(8c) 



Equations © define the classical geometric pic- 
ture associated to quantum two-level systems: the 
unit vector S precesses around the vector B just 
like a classical gyromagnet precesses in a magnetic 
field. By means of a canonical transformation, S = 
( \Jl — q 2 cos p, yl I 2 sm P,—q), we get the following 
bi-dimensional canonical representation 



H = - (B x cosp + By ship) \J\ - q 2 - B z q, (9a) 
q = (B x sinp- B y cosp) \j\ - q 2 , 



p = — (B x cosp + By sinp) 



(9b) 
- B z , (9c) 



in which the amplitudes B x y ^ z arc functions of time. The 
square roots have as consequence that © does not sat- 
isfy a Lipschitz condition (see, e.g., Q, Theorem 1.2) and 
thus there is not, a priori, a local existence and unique- 
ness theorem (see, e.g., Q, Theorem 1.1). Thus, a com- 
parison with the (nonsingular) dynamics l|8bjl is neces- 
sary to select the correct solutions. Representation 
is, however, important for the construction of the stro- 
boscope maps, which are numerical tools to attest the 
integrability of the system. 
We have now: 

Result 2 : The system (8) is intcgrablc. 

In order to show this we consider the classical distance 

function on the Bloch sphere 



£>(i) = ||Si-S 2 



(10a) 



where different indices refer to different initial conditions. 
By 0, it follows from (a-E)(b-E) = (a- b) + tT, ■ (a x b) 
that 



D(t) = ^2Tv( Pl -p 2 ) 2 , 



(10b) 



where the sub-indices refer to different initial con- 
ditions on the Schrodingcr wave-function pi, 2 {t) = 

101,2(^)^1,2^)1, With 



Pi,2(t) = U(t) Pl , 2 (0) U\t), 



(11) 



where U is the unitary propagator satisfying iU = 
H{t)U(t). Equation l|10a|) allows us to define a Lyapunov 
exponent as 



A = lim lim — In 

D(0)^0 t-»t» t 



D(t) 



D(0) 



(12) 



However, by l|10b|) and i|ll|) it is immediate that D(t) = 
D(0) and thus, by (JT2J), A = 0. When a full set (i.e., 
equal to the number of degrees of freedom) of constants 
of the motion in involution does not exist, so that we 
have a system which is not integrable, sensivity to ini- 
tial conditions is certain to exist at least in some region 
of phase space, i.e., there is inevitably chaotic behavior. 
Thus we may regard the proof that A = as a proof of 
integrability, although it may be difficult to find explic- 
itly an additional constant of the motion for the general 
autonomous system equivalent to © (see, e.g., [llj, for 
the description of these autonomous systems) . It is worth 
emphasizing that our proof of integrability is valid for an 
arbitrary time-dependent field, including the quasiperi- 
odic case. This contradicts the basic assertion of ^3 
and confirms the result of jl3j. More importantly, we 
show the real reason for the integrability observed in 01 
using double Poincare sections. 

Equations (JSJ describe the interaction of one two-level 
atom with the quantized radiation field (q.r.f.) in the Flo- 
quet approximation, where the average number of pho- 
tons becomes very large, keeping the photon density con- 
stant . It may also be regarded as a special case of the 
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(mean field or classical) limit of the expectation-values 
of the Heiscnberg equations of motion for a system of 
many atoms interacting with the q.r.f. in suitable (e.g. 
coherent) states 01 • I n this limiting process, Result 2 
is no longer generally true and, indeed, in the single- 
mode case, in which the electric field is periodic in time, 
the resulting equations do exhibit chaotic behavior for 
the model with an initial population, when the polariza- 
tion (i.e. the macroscopic transition dipole moment) is 
nonzero [lfij . Only the special situation of zero polariza- 
tion p reduces to our equations |JSJ, in which case it is 
wrongly asserted in that the model is integrablc due 
to condition JHcJ). Indeed, l|%c|) solely reduces the three 
spin variables to two independent ones, i.e., one obtains 
a one-degree of freedom system in an external periodic 
field. The equivalent autonomous system is thus in gen- 
eral chaotic due to the nonlinearity in the time-variable 
[i"l"|| . It is our Result 2 which shows that integrability 
holds as a consequence of the unitarity of the underlying 
quantum dynamics. Finally, the full one-mode quantum 
Hamiltonian has striking similarity to our case, because 
it is unitarily equivalent to a Hamiltonian where the q.r.f. 
and the atom are decoupled |l7|. and thus is integrablc 
in this sense. The reason of this similarity is that the 
Floquet approximation preserves some of the quantum 
structure |14j . while the classical approximation (limit) 
is rather singular [l^. Only in the case p = do both 
approximations agree. 

In the special case (|5J it is easy to prove integrability 
explicitly. Consider the following generating function and 
respective canonical transformation: 



F 3 {p,Q,t) 



P 



-pQ + utQ, 
op 



dQ 



dF* 



/C(Q,P) = H + -^=H + luQ. 



(13a) 
(13b) 

(13c) 

(13d) 



The new Hamiltonian K, is autonomous and the sys- 
tem is integrablc. Moreover, a stroboscopic map in 
the original space (p, q) can be constructed at strobo- 
scopic instants t k = 2Trk/tu, with k G N, by noting that 
(PkiQk) = (Pkilk)- Since K is a constant of motion, 
its contour curves will correspond to the stroboscopic 
map for irrational values of B/uj, being B = 2\JBq 2 + £1 2 
the characteristic frequency of the integrated motion and 
fl = Bi - lu/2. By (0), HSU and JT3J) we get the contours 
for the R stroboscopic map, namely 



K = 2B ^/l 



q k z cosp k 



2(B 1 --)q k . 



(14) 



From now on we use the notation r(tfe) = ?*fc- A further 
interesting result for this system is the analytical corre- 
spondence between periodic orbits in the classical phase 
space and eigenstates in the Hilbert space . 



For the NR system , the absence of a rotation sym- 
metry destroys the manifest integrability of the model. 
Notice that the modulus of the magnetic field is no 
longer constant and finding a rotating frame does not 
suffice. Numerical calculation is now required. For 
convenience, hereafter we shall work with Unr = 
— 2(i?3, 0, — i?2 cosojt), which results from a rotation of 
the field through the angle 7r/2 around the y-axis. 
We also set B3 — > Bq and B2 — > —B\ in order to es- 
tablish contact with the field (J2J). A first remarkable 
result found numerically is the linear relation between 
the time-dependent energy, H(t). and the coordinate q(t) 
at the stroboscopic instants [l8j. It allows us to write 
Hk = £(9o,Po) ~7Qk and, consequently, 



£(go,Po) = \/l - q k 2 cosp k - 2 [Bi 



7 



fffc, (15) 



which is the equation for the NR contours. Comparison 
between l|14|) and i|15|) emphasizes the resemblance be- 
tween NR and R dynamics at stroboscopic instants. It 
illustrates the integrability of the NR system. The qual- 
itative agreement between NR and R maps are shown in 
Fig^ which was made as follows: we determined 7 nu- 
merically by fittings in the NR system and then we used 
a R field with uj = 7. 




FIG. 1: (a) Stroboscopic map for the NR system with Bo = 
1.0, Bi = 1.5 and u — 3.0 and (b) the corresponding map for 
the the R system with w = 7 = 4.9559 (see eq. 1141 . This 
value was obtained by the fitting procedure [13 . 

The contour pattern l|15|) "chosen" by the NR dynam- 
ics is not obvious. In fact it cannot be inferred directly 
from the Hamiltonian for the NR field. The choice of 
one specific contour in the set of possibilities offered by 
the time-dependence of the system suggests the existence 
of an average dynamics which underlies integrability (see 
[3 for details). 

The integrability of the system is made evident via 
averaging theorem in other two limiting regimes, namely, 
the strong coupling (Bo <C u) and the resonant weak 
coupling (2Bq = u) and B\ <C u>) regimes. The validity 
of the RWA is also numerically seen to extend beyond 
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the time scale 0(1/ B^) 0. Thus, integrability and its 
consequences explain the wide validity of the RWA (see 
also |l9j for numerical results). 

Last but not least, we describe our results for the NOT 
operation in the NR system. Although analytical solu- 
tions were not found in this case, we may use (|14|) and 
(|15(l to make predictions. 

General principles of quantum computation require 
that the NOT operation be unitary |2flj. References 
HH E3 propose different constructions of the universal 
unitary NOT operation, which takes as input a qubit in 
an arbitrary state \ip) and generates an output as close 
as possible to its orthogonal state \ipx). Since we are 
primarily concerned with the dynamics, we consider the 
special cases where a unitary operation leading from \ip) 
to \i/j±) exists. 

Following the derivation of (|10bjl . one may obtain 
2\(ip \i; t )\ 2 = 1 + S • S t . The NOT operation occurs 
when the dynamics produces, at the instant t not , a state 
\ipnot) orthogonal to the initial state \ipo)i such that 
(V'ol^/'noi) = 0. Classically, it corresponds to S not = — Sq. 

According to J7J), the corresponding initial quantum 
states which allow the NOT operation are simply ob- 
tained by putting (po, qq) given above into 



1 



-e ip \ 



(16) 



which emphasizes that the canonical momentum p plays 
an important role of a relative quantum phase. In |l8| we 
have found analytically the regimes for the realization of 
the unitary NOT operation on the R system. In particu- 
lar, when the resonance condition uj 2 = Bq 2 + (Bi—uj/2) 2 
is satisfied, the NOT operation will occur periodically at 
tnot = (2n + l)ir / to for initial conditions given by arbi- 
trary po and qo = 0. Very interestingly, this regime also 
occurs for the NR system. Following the analogy between 
7 and u> given by (|14|) and 1|15|) we determine numerically 
the parameters satisfying the required resonance and the 
universality class of initial states for the NOT operation 
in the NR system (see Fig[2J). 

The predictions of NOT operations in the NR sys- 
tem were based on comparisons with stroboscopic maps 
of the R system. These classical mathematical tools, 
based on the concept of trajectories flowing in the phase 
space, were crucial in the determination of the connec- 
tion between the R frequency u> and its NR counterpart 
7. It stresses the usefulness of the classical gyromagnet 
picture, being an example of classical analysis allowing 
quantum predictions. 

R.M.A. thanks financial support from FAPESP 
(Fundacao de Amparo a Pesquisa do Estado de Sao 
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FIG. 2: NR NOT operation for w = 1.0, 7 = 1.486, B 1 = 1.5 
and Bo = 1.279 occurring at t no t = 57r for several initial 
conditions given by (po,qo) = (f,V), which corresponds to 
a universality class of quantum states with purely complex 
relative phase. 
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